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ABSTRACT 


This report describes a series of low-speed airfoil designs based 
on modification to the NACA 64-206 and 64^-212 airfoils. Designs are 
based on potential flow theory. The report describes one of a series 
of airfoil modifications carried out under Contract NAS 2-8599, Appli- 
cation of Multivariable Search Techniques to Optimal Wing Design in 
Non-Linear Flow Fields. Mr. Raymond Hicks of National Aeronautics and 
Space Administration’s Aeronautical Division, Ames Research Center, served 
as contract monitor for the study. 
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THEORETICAL EFFECT OF MODIFICATIONS TO THE UPPER SURFACE OF 

TOO NACA AIRFOILS USING SMOOTH POLYNOMIAL ADDITIONAL THICKNESS 
DISTRIBUTIONS WHICH EMPHASIZE LEADING EDGE PROFILE AND 
WHICH VAR V LINEARLY AT THE TRAILING EDGE 

by Donald S. Hague and Antony W. Merz 
Aerophysics Research Corporation 

SUMMARY 

An investigation has been conducted on the Lawrence Radiation Center, 
Berkeley, CDC 7600 digital computer to determine the effects of additional 
thickness distributions to the upper surface of the NACA 64-206 and 
64^-212 airfoils. The additional thickness distribution had the form of 
a continuous mathematical function which disappears at both the leading 
edge and the trailing edge. The function behaves as a polynomial of 
order at the leading edge, and a polynomial of order z^ at the trailing 
edge. In the present study, is a constant and is varied over a 
range of practical interest. The magnitude of the additional thickness, 
y, is a second input parameter, and the effect of varying and y on 
the aerodynamic performance of the airfoil was investigated. Results 
were obtained at a Mach number of 0.2 with an angle-of-attack of 6° on 
the basic airfoils. All calculations employ the full potential flow 
equations for two dimensional flow. The relaxation method of Jameson is 
employed for solution of the potential flow equations. 

As was found in earlier investigations of these airfoils, using other 
types of additional thickness distribution, increasing the additional 
thickness, y, tends to increase both the lift and the adverse pitching 
moment cofficients. In the present investigation variations in the shape 
parameter, e^, also change the lift and moment coefficients proportionally 
and monotonically. For the range of coefficients examined, the lift 
coefficient was nearly insensitive to changes in e , while the pitching 
moment coefficient showed stronger variations to this parameter. For 
both the 64-206 and 64^-212 airfoils, however, the lift and moment 
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coefficients were much more sensitive to changes in y than to changes in 

V 

Additional thickness can be made to produce significant reductions 
in peak pressure coefficient. This is particularly true for the 64-206 
airfoil, which in its unmodified form has a very high pressure peak at 
the leading edge, due to the small radius of curvature at this point. 

For this airfoil, the present results indicate a complex dependence of 
the peak pressure coefficient on the parameters and y. Study of this 
dependence was limited due to the increases in adverse pitching moment 
which accompanied the decreases in peak pressure coefficient. For the 
limited range of parajqpter variations conducted, the peak pressure and 
lift coefficient of the 64^212 airfoil followed similar patterns with 
variations in e 1 and y. For this airfoil the peak pressure reduced 
monotonically with e^, while the lift increased monotonically with 
thickness, y. 

It should be noted that viscous effects are neglected in the present 
analysis. At the higher lift coefficients the effect of viscosity could 
be significant. Further investigations incorporating a viscous flow model 
are therefore desirable. 


INTRODUCTION 

The National Aeronautics and Space Administration and others are 
currently conducting a series of theoretical and experimental studies 
to define airfoil sections having improved performance from the aspects 
of lift, drag, pitching moment, or pressure distribution characteristics, 
references 1 and 2. Analytic investigations using airfoil surface repre- 
sentations based on high-order polynomials may result in impractical 
profiles; for example, very thin trailing edge thickness distributions 
or severe reflexes in the profile. The present study employs a continuous 
polynomial arc having two free parameters, whose characteristics are 
selected to avoid such problems. Previous optimization studies using mult 
variable search techniques, references 1, 4 and 5, generally indicate that 
shape changes which provide increased lift produce unfavorable changes in 
moment characteristics. Conversely, profile changes which improve the 
moment characteristics decrease the lift coefficient. With the low-order 
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model of the present investigation a systematic examination on the effect 
of profile changes can be carried out. This was accomplished and the 
trends previously revealed by optimization studies were confirmed. 
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MATHEMATICAL MODELS 

Potential Flow Equation 

Potential flow analysis is based on solution of the two-dimensional 
potential flow equation 

(a 2 -u 2 ) ^ + (a 2 -v 2 ) <}> yy - 2uv «j> xy = 0 
where <t> is the velocity potential, u and v are the velocity components 

U = d> , V = A 

' x* v y 

and a is the local speed of sound determined from the energy equation and 
the stagnation speed of sound 

2 2 ,v - U / 2 2. 

a = a 0 - C- T - — 0 (u + v ) 

Solutions are obtained by Jameson's finite difference scheme, reference 6. 


AIRFOIL PROFILE REPRESENTATION 


Basic Airfoil 
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Ordinates for the basic NACA 64-206 and 64^-212 airfoils were approximated 
by four cubic chain polynomials in the manner of Hicks 

2 3 

y. - a F. + a. x + a x + a x ; j = 1,2, 3, 4 
J 1 J 3 j J 

Coefficients in the four polynomial arcs are selected on the following 
basis : 

j = 1 - Arc represents forward portion of upper surface 

F i - 

. j=2 - Arc represents aft portion of upper surface 

f 2 = 1 

j=3 - Arc represents forward portion of lower surface 


F,= V' 
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j o 4 - Arc represents aft portion 

V 1 

The coefficients a. are determined by introducing four boundary conditions 
x j 

on the airfoil profile in each of the four airfoil arcs, Crout's method 
for triangularization and back substitution is used to solve the resulting 
system of linear equations. Note that if four points are specified on the 
aft portion (i = 2 or 4) , a discontinuity in slope occurs where the poly- 
nomials join. This produces a small ripple in the pressure distribution 
at the juncture point. However, since the juncture occurs at a region 
of small slope (x = .$} the effect is not significant. 

Computer-generated plots of the NACA 64-206 and 64^-212 airfoils, 
together with the associated pressure distributions predicted by potential 
flow theory, are shown in Figure 1. Table I lists the polynomial coefficients, 

a , which are used for the basic airfoil representations. 
i . 


Additional Thickness 


The upper surface of the basic airfoil is modified by addition of 

the thickness-distribution function, 

e l e 2 

Ay(x) = Ax (1 - x) 

where = 1, for the present study, It is shown in Appendix A that the 
magnitude parameter. A, can be expressed in terms of the maximum thickness, 
y, by the equation. 


A = 


y 


Ci + El ) 



Representative functions Ay(x) are shown in Figure 2, for e x = .25, .50 and 
.75. Notice that the slope of the additional thickness distribution is 
infinite at the leading edge (x = 0) and is finite and positive at the 
trailing edge (x = 1) . Notice also that the leading edge shape is very 
sensitive to this exponent, and that values of e.^ less than .25 produce 
a very blunt nose on the airfoil. The chordwise location of the point of 
maximum thickness is shown in Figure 2(d) as a function of the shaping 
parameter e.. 
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TABLE I. COEFFICIENTS FOR AIRFOIL REPRESENTATIONS 


64-206 Airfoil 


Upper Surface, Forward 
Upper Surface, Aft 
Lower Surface, Forward 
Lower Surface, Aft 


a 0 

a l 

a 2 

a 3 

.07784 

.00098 

-.01209 

-.10986 

.02453 

.10936 

-.18959 

.05569 

-.06050 

.05728 

.08033 

,14000 

-.01026 

-.09394 

.21987 

-.11567 


64^212 Airfoil 


Upper Surface, Forward .14641 -.01941 .00039 
Upper Surface, Aft .03227 .27010 -.50591 
Lower Surface, Forward -.12650 .07374 -.09721 
Lower Surface, Aft -.01021 -.29552 .60784 






OPTIMIZATION STUDIES 


In previous airfoil optimization studies, (references 4 and S), the 
modifications to the upper surface took the form of a pair of quadratic 
arcs, which were cotangential at the point x, y. These parameters are 
respectively the chordwise location of the maximum addition thickness 
and the value of this thickness. Both lift coefficient and moment 
coefficient were considered as performance indices in this development. 
The present study is also concerned with a two variable optimization 
problem using the leading edge thickness distribution exponent, e^, and 
the magnitude of additional thickness, y. 

Lift Coefficient Maximization 


In general, maximi.gav.ion of lift coefficient has the form 


where 


<P = Max 


= i)Ap(x)dx 


and the integration is around the airfoil contour. The airfoil contour 
in the present study and those of references 4 and 5 are completely described 
in terms of two parameters, and c^* For the present airfoils 


<P = Max [cj = Max C Jj (E 1 ,y)J h Max 


where 


This two variable optimization problem can be solved by use of multi- 
variable search techniques, for example, a combination of directed random- 
ray and pattern searches, references 3 and 7. 

Examples illustrating this type of search procedure have previously 
been presented in references 4 and 5. However, the low dimensionality of 
the present problem (two parameters) permits the solution of optimization 
problems by inspection of graphical results. This procedure is employed 
in the present report. Other optimization problems of interest are described 
below. 




Moment Coefficient Minimization 


Minimization of the moment coefficient has form 


wnere 


♦ “ Minj^j “ Min [ C M Cg 1'^] 5 Min [VV 0 ^] 


C M M(x - 1/4) Ap(x)dx 


In previous studies moment minimization resulted in a solution directly 
opposed to lift maximization. The position of maximum thickness moved 
forward and the amount of additional thickness was minimized. Thus in those 
studies the basic airfoil had less adverse moment than any airfoil generated 
by addition of the specified thickness to the upper surface of the airfoil. 


Other Optimization Criteria 

Other airfoil performance criteria can be considered, which typically 
involve compromises between lift, moment and peak pressure coefficients. 
Such modified criteria can take any of the following forms: 


1. Maximize a linear combination of lift and moment: 


= Max 

2. Minimize the moment at a specific value of lift: 


* 

v a y 


* 


= Min C u 

M _ 


3. Minimize the peak pressure at a specific value of lift: 


* 


Min 



In the studies reported in references 4 and 5, the parameters avail- 
able for airfoil modification (x and y) , both were varied over a large 
range. This permitted a straightforward interpretation of results, such 
that optimal parameter pairs for a given performance criterion could be 
determined by inspection. As noted above, this procedure is also followed 
in the present study. Free variables for the present study are the 
parameters and y. 


SYSTEMATIC AIRFOIL SHAPING 


The present study is primarily concerned with a limited but systematic 
investigation on the effects of varying the leading edge thickness magnitude 
and distribution for the NACA 64-206 and 64^-212 airfoils. The parameters 
y and e^ are varied over ranges which are sufficient to permit qi tative 
conclusions as to their effects on lift, moment, and peak pressure coefficients. 
For each pair of such parameters, a plot of the airfoil and its associated 
calculated pressure distribution are given in Figures 3(a) to 3(r). 

The first 12 of these plots relate to modifications of the 64-206 airfoil, 
and the last 6 are related to modifications of the 64^206 airfoil. These 
pressure signatures differ from those of the basic airfoils (Figure 1) 
chiefly in the magnitude of the peak pressure at the leading edge. Increases 
in both Ej and in y tend to soften the pressure variations over the upper 
surifp by reducing the leading edge peak and by increasing pressures over 
t,hs? central and trailing edge regions. In the case of the modified 64-206 
airfoil high values of ultimately reverse this trend and the strong 
overpressure peak reappears. 

Lift and moment coefficient variation for the two airfoils are shown 
in Figure 4, and the effects of varying the parameter and y are apparent. 

For both airfoils, the lift increases only slightly with the exponent 
while the thickness y has a more pronounced influence on the lift increment. 

The adverse moment coefficient also rises sharply with additional thickness, 
while the exponent has a more significant influence on C^. Combined 
lift vs. moment results are given in Figure 5, which shows that the least 
adverse moment is obt?iu'-u %Z ? given y, or at a given C^, with the 
smallest value of e^. This to*"-: londs to a relatively blunt leading edge 
on the airfoil. 

Variations of the peak pressure with the parameters e^ and y are shown 
in Figure 6 for the two ?,irfoils being studied. The pressure variation 
for both airfoils is minimized at a particular y by a specific choice of 

For the 64^-212 airfoil peak overpressure at a given lift coefficient 

is minimized by using the largest leading edge exponent value, e^. The 

Ci-C variation is more complex for the 64-206 airfoil, however, in 
L p 
*max 

that the constant e 1 loci cross each other. The envelope of these curves 
is given in Figure 7, and it defines the minimum peak pressure for a 
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given lift. The small number of data points available in this study does 
not permit accurate cross-plotting, so the estimated minimum pressure peak 
values are shown in the cross-hatched area. Despite the uncertainties in 
the cross-plotting procedure it is evident that the present family of 
airfoil designs produce lower peak overpressures for a given than the 
airfoils previously obtained through biquadratic modifications. 

Qualitative results of this parametric study are summarized in Table II. 
These conclusions follow directly from the results shown in Figures 4 to 7. 
Figure 7 also presents the minimum peak overpressures obtained with biquadratic 
airfoil modifications in references 4 and 5. 

CONCLUSION 

A numerical study has been completed of a class of modifications to the 
NACA 64-206 and 64^212 airfoils. Systematic changes in the upper surfaces 
of these airfoils were studied by independent variations in the thickness 
und leading edge thickness distribution exponent. The Mach number and 
angle-of-attack were constant during the study, and the results are summarized 
as follows: 

1. Pressure distribution is moderately sensitive to leading edge 
profile and to additional thickness. 

2. Lift coefficient is nearly independent of the leading edge profile, 
but increases with additional thickness. 

3. Adverse pitching moment increases with additional thickness and 
with increases in the leading edge profile exponent. All 
modifications to the airfoils increased this adverse moment. 

4. Peak pressure for a given lift coefficient can be considerably 
reduced by careful selection of the leading edge additional 
thickness distribution exponent. Somewhat lower peak pressures 

at a given lift are possible using the present airfoil modifications, 
as compared with the "biquadratic” modifications of references 4 
and 5. 

5. For a given lift coefficient, adverse pitching moment is minimized 
by reducing the leading edge profile exponent. 


Criterion 


Exponent Cj 


Thickness y 
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Max 


Min 


64*206 
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Max 


Min 


Max 


Min 


Basic airfoil is best 















AND PRESSURE DISTRIBUTIONS 
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FIGURE 3(a). MODIFIED 64-206 AIRFOIL y - .03, Cj - .10 
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FIGURE 3(b). MODIFIED 64-206 AIRFOIL y - .03, e. - .15 













FIGURE 3(0). MODIFIED 64^212 AIRFOIL y - .03, Cj ■ .5, c 2 ■ 1 



FIGURE 3 (p) . MODIFIED 64^212 y » .06, Cj - .2, c 2 » 1 
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(b) 64^-212 Airfoil 


FIGURE 5. LIFT AND MOMENT VARIATIONS 
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FIGURE 6. LIFT AND PEAK PRESSURE VARIATIONS 




APPENDIX A 
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CHORDWISE LOCATION OF MAXIMUM THICKNESS 


The distribution function used in this study is 

C 1 c 2 
Ay(x) * Ax (1 - x) 

and the variation of this function is smooth for 0 s xs. 1. The 
point of maximum additional thickness occurs when 


(A-l) 


Ay' (x) - Ay(x) 


-1 

CjX 


c 2 (1 - x) 


>] .0 


This shows that the chordwise location of the point of maximum 
thickness is 


x * 


1 


€ 1 * e 2 


and the value of the maximum thickness is then found in terms of 
the parameters as 


Ay, 


max 


r C 1 e 2 
C 1 C 2 


L (e l * e 2 5 


e l * e 2 


(A- 2) 


(A-3) 


( A- 4 ) 


For the case studied in this report, e 2 * 1. and the parameter 
A is 

1 ♦ e, 

C i * 

e l 


1 


(A-5) 
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